SCIENTIFIC 

REPpRTS 



»H5 




SUBJECT AREAS: 

STRUCTURE OF SOLIDS 
AND LIQUIDS 

METALS AND ALLOYS 

MECHANICAL PROPERTIES 

MINERALOGY 



Microstructural fingerprints of phase 
transitions in shock-loaded iron 

S. J. Wang', M. L Sui 2 , Y. T. Chen 3 , Q. H. Lu\ E. Ma 4 , X. Y. Pei 3 , Q. Z. Li 3 & H. B. Hu 3 

1 Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 
1 1 001 6, China, institute of Microstructureand Property of Advanced Materials, Beijing University of Technology, Beijing 1 001 24, 
China, institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 62 1 900, China, 4 Department of Materials 
Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA. 



Received 
4 October 2012 

Accepted 
11 December 201 2 

Published 
18 January 201 3 



Correspondence and 
requests for materials 
should be addressed to 
M.L.S. (mlsui@bjut. 

edu.cn) 



The complex structural transformation in crystals under static pressure or shock loading has been a subject 
of long-standing interest to materials scientists and physicists. The polymorphic transformation is of 
particular importance for iron (Fe), due to its technological and sociological significance in the development 
of human civilization, as well as its prominent presence in the earth's core. The martensitic transformation 
a-^z (bcc^hcp) in iron under shock-loading, due to its reversible and transient nature, requires non-trivial 
detective work to uncover its occurrence. Here we reveal refined microstructural fingerprints, needle-like 
colonies and three sets of{112}<lll> twins with a threefold symmetry, with tell-tale features that are 
indicative of two sequential martensitic transformations in the reversible a-^z phase transition, even 
though no £ is retained in the post-shock samples. The signature orientation relationships are consistent 
with previously- proposed transformation mechanisms, and the unique microstructural fingerprints enable 
a quantitative assessment of the volume fraction transformed. 

At ambient conditions, Fe is stable in its oc-phase with the body-centered-cubic (bcc) crystal structure. At 
high pressures, Fe takes on the 8 form, switching to the hexagonal-closed-packed (hep) lattice. While the 
phase diagram of Fe under hydrostatic pressure is well established 1 , it is much more challenging to 
ascertain what happens in Fe when it is subjected to shock-loading. Understanding dynamic phase evolution 
is critically relevant to many applications of iron and steels under explosion and shock processing conditions. 

Nearly sixty years ago, an oc^s martensitic transformation (MT) was postulated for Fe upon shock-loading (to 
13 GPa) 2 ' 3 . But the confirmation of this shock- induced transformation requires non-trivial detective work, 
because the high-pressure s-phase induced during the loading transient transforms back into the a-phase 
immediately after the shock is over. As a result, the reversible oc^s transformation leaves no s-phase behind, 
for direct observation in the microstructure of postmortem samples. Earlier experiments had to rely on indirect 
evidence such as wave profile measurements 2 " 7 : the stable three-wave shock fronts of the wave profiles measured 
during shock-loading suggested a polymorphic phase transition. More recent efforts succeeded in acquiring 
diffraction patterns from the s-phase via in-situ x-ray diffraction both in single crystal and polycrystal 8 " 10 , and 
captured signatures from in-situ extended x-ray absorption fine structure 11 in laser shocked Fe. These experi- 
mental findings are consistent with the predicted a^s transformation in recent multi-million-atom molecular- 
dynamics simulations 912 " 15 . However, key questions remain open in this area: 1) While in-situ X-ray diffraction 
measurements provided information about the atomic pathway of the phase transition, and the width of diffrac- 
tion peaks could be used to deduce the microstructural scale involved, the previous experiments, including those 
that examined recovered samples after shock 1617 , were not able to directly image the morphological features (such 
as size, shape and density of domains, defects, twins, interphase boundaries, etc.) that are characteristic or 
indicative of the oc-^s transition. 2) Without such microstructural information, the experiments so far cannot 
quantitatively assess the volume fraction actually transformed to the s-phase, in different locations of the shocked 
sample, especially when the loading geometry/condition is complex and the sample is of bulk form. In recovered 
samples, only a minute amount of s-phase could be retained by rapid quenching using a femto-second laser 18 . 

The present work demonstrates that insight into these issues can be obtained by performing ex-situ transmis- 
sion electron microscopy (TEM) experiments, on recovered samples of shock-compressed Fe. In particular, we 
have detected under TEM signature microstructures retained in the post- shock sample, including needle-like 
colonies arising from shock- induced a^s transition and three sets of ultrafine twinned a domains inside the 
colonies which originate from the three-fold prismatic planes of the s phase during the s^a reverse transition 
upon unloading. These microstructural fingerprints can be used to infer that the a^s^a MTs have indeed 
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Figure 1 | Inverse pole figure color map overlapped with twin boundaries in a-Fe. (a) before and (b,c,d) after shock-loading. Many twin lamellae appear 
after shock-loading. The {112}<111> ( J]3) twins with boundaries delineated by black lines are dominant, while the {332} < 1 13> twins (white lines) are 
very limited. The colors in the IPF map correspond to the crystallographic orientation normal to the observed plane, indicated in the stereographic 
triangle superimposed in (a). The insets in (b,c) are the enlarged images of the corresponding regions marked by arrows. 



occurred in the shocked material, with unique features that are in line 
with specific mechanisms proposed for both the forward and reverse 
a-s MTs. More importantly, based on the postmortem microstruc- 
tures, we can assess quantitatively the volume percentage trans- 
formed under different shock-loading conditions and in different 
regions of the sample volume. This investigation thus provides 
quantitative information that is essential to the development and 
validation of physical models and computational algorithms (such 
as macroscopic treatment via numerical calculations) about shock- 
induced phase transformations. 



Results 

Bulk polycrystalline Fe with a hollow cylindrical shape was shocked 
by explosive loading on outside surface (see the details of experi- 
mental setup in the Methods section). The peak pressure near the 
loading surface has been estimated as —16.4 GPa (higher than the 
pressure threshold anticipated for the a-^s MT 2 ' 3 ). The original 
sample consists of equi-axed a-Fe grains with sizes in the range of 
tens of micrometers, as indicated by the electron backscatter diffrac- 
tion (EBSD) orientation map in Fig. la. Fig. 1 shows inverse pole 
figure (IPF) color map overlapped with twin boundaries (black and 
white lines). The colors in the IPF map correspond to the crystal- 
lographic orientation normal to the observed plane, indicated in the 
stereographic triangle superimposed in Fig. la. After shock-loading 
from the cylinder outside surface, fine sub-grain features can be 



observed in the original a-Fe grains, as seen in the IPF color maps 
in Fig. 1, as well as in the TEM micrographs in Fig. 2. In what follows, 
these microstructural details will be shown to be revealing regarding 
phase transitions, while X-ray diffraction merely confirmed that a-Fe 
remained the only phase present after the shock throughout the 
sample, with no indication of phase transformation. 

Specifically, Figs, l(b-d), taken at locations about 4.6 mm (local 
pressure is calculated as ~ 12.7 GPa by hydrodynamic code shown in 
the Methods section), 1.9 mm (-13.3 GPa) and 1.2 mm (-14.2 GPa) 
inward from the loading surface, respectively, show ultrafine lamellar 
structures. Based on the orientation information revealed by EBSD 
orientation map, these planar defects are twins and their boundaries; 
see TEM details later. There are more and finer twin lamellae near 
the cylinder outside surface [Figs. l(c,d)] than that much further 
away from the surface (Fig. lb). 

The twinned structures near the outside surface have been inves- 
tigated carefully under TEM. Fig. 2a is a montage collection of 18 
bright-field TEM images, taken one by one from an area of about 
35X35 |im 2 (extracted from one large original a-Fe grain in 
Supplementary Fig. SI) viewed along the <110> zone axis. In this 
large- area TEM view, in addition to some dislocation/strain contrast, 
there are two levels of salient microstructures that will be discussed in 
detail below. On the finest level, there are numerous fine lamellar 
structures and boundaries, which are the twin lamella features men- 
tioned earlier (see EBSD findings in Fig. 1). A close-up view is shown 
in the high-magnification TEM micrograph in Fig. 2b, which 




Figure 2 | Microstructural fingerprints of phase transitions in shock-loaded iron, (a) TEM micrograph showing the microstructure of shock- 
compressed Fe. (b) Enlarged view of the yellow rectangle area in a): TB^ TB 2 and TB 3 (dashed, dotted, and dash-dot green lines) stands for three sets of 
{112} twins; T 332 and T( 332 )' (red solid lines) for two sets of {332} twins. The triangle with black solid lines shows the included angle of the three sets of 
{ 1 12} bcc twins, (c) SAED of the three sets of { 1 12} twins, (d) SAED of T 332 . (e) Contours of the boundaries of needle-type-martensite-like regions in a) 
(white dashed lines). 
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enlarges the yellow boxed region in Fig. 2a. One observes three sets of 
straight edge-on boundaries (marked with green lines) of intersect- 
ing lamellar bands. These, as seen from the selected area electron 
diffraction (SAED) patterns in Fig. 2c, are three sets of {112}<111> 
twins with a threefold symmetry. In addition, Fig. 2d gives the SAED 
pattern indicating the {332} < 1 13> twin relationship for the bound- 
aries (red solid lines) at the upper right corner of Fig. 2b. These 
{332}<113> twins coexist with the dominant {112}<111> twins. 
But since they also exit in plastically deformed regions far away from 
the loading surface (Fig. lb), they cannot be taken as a signature 
indicative of phase transformation 17 . 

At a coarser level, a salient microstructure in Fig. 2a is the presence 
of many needle-like regions (for the clear boundaries between col- 
onies and matrix, see the original TEM image in Fig. S2), enclosed by 
white dashed lines. The needle-like colonies form networks, dividing 
the original single a-grain into many blocks. A cartoon delineating 
the network is displayed in Fig. 2e. The multiple sets of twin lamellae 
discussed above are contained inside these needle-like regions. 

To recapitulate, both before and after shock-loading the Fe is in its 
a-phase. But a hierarchical microstructure composed of two levels/ 
types of unusual features has been detected after the shock. These 
salient features are the fingerprints that memorize and reflect the 
(phase transformation) history the material has been through. In 
the following, we explain the origins of these features and associate 
them with the phase transformations experienced. 

We first discuss the edge-on boundaries (marked with three types 
of green lines in Fig. 2b). As shown above, these are three sets of 
lamellae of {112} twins, in three different orientations. As seen in the 
micrograph, they form an acute angle of 60° in between when they 
intersect (inset of Fig. 2b). To explain this and relate it to the a-^s-^a 
MTs, we first point out that when viewed along the < 1 10> zone axis 
of the bcc a-Fe, two {112} planes should have formed an inclusion 
angle of 70.5° or 109.5°, rather than 60° or 120°. Therefore, the three 
sets of 60° {112} planes in Fig. 2b are not {112} twins of a-Fe within a 
single crystal (grain). In the meantime, we note that martensite 
structures are well known to often comprise a twinned microstruc- 
ture 19 , as it helps to relieve the strains that accompany the MT. The 
three sets of twinned lamellae are thus three variants of a martensite, 
forming three domains that later impinge and intersect. In other 
words, each set of the {112} lamellar twins arises and grows out of 
a separate domain (subgrain); and when two neighboring sets 
impinge, they meet at the 60° angle. 

What, then, is the origin of the 60° orientational relationship and 
threefold symmetry (Fig. 2b)? The answer is that the martensite 
formation is not from the bcc a-phase but rather from an intermedi- 
ate phase, which has an intrinsic three-fold symmetry. In other 
words, before the MT mediated by twinned structures, the region 
in question (the needle-like colony) should have already undergone 
the a-^s transformation. In the hep s-phase, MT habit plane is 
expected to be {l-100} hcp which has six-fold symmetry 20 . There are 
three crystallographically equivalent {1-100} planes, and this three- 
fold symmetry is retained in the reverse transformation back to the 
a-phase. 

To lend support to this assertion, we note that Zhang et al 20 have 
observed earlier a similar microstructural feature for a-Ti(hcp)^P- 
Ti(bcc) MT induced by rapid heating: three sets of edge-on {112} 
martensite twin lamellae of P-Ti(bcc) with an angle of 120° (or 60°), 
surrounding an untransformed a-Ti(hcp) region in a triangle shape. 
The orientation relationship between a-Ti and P-Ti is determined to 
be (l-100) a . Ti //(-112)p. Ti (and [0001] a _ Ti //[110] p _ Ti ). Similar to 
these previous findings, the three sets of {1 12} ot _ Fe lamellar twins 
in Fig. 2b should be the a-Fe martensite variants (lamellar twins) 
originated from three different {l-100} 8 . Fe planes, and the angle of 
60° is inherited from the three-fold symmetry characteristic of the 
{l-100} 8 _ Fe planes. Each of the three variants grows out of one pris- 
matic {1-100} plane of the s-phase, thus meeting with 60° angles. 



The microstructure observed is therefore a signature revealing the 
a-^s^a transformation sequence. Shock-loading does induce the 
a-^s transition, and the intermediate s-phase transforms back to a- 
phase during the ensuing unloading. A similar argument (see the 
Discussion section) can easily exclude the involvement of an inter- 
mediate face-centered-cubic y-phase 13 ' 21 in the MTs. 

We next discuss the networked needle-like regions (outlined by 
white dashed lines in Fig. 2a). If the bundles of the twin lamellae 
contained inside the needles come from the s^a reverse transition 
as discussed above, then the needle-type regions themselves would be 
the locations that have experienced the a^s forward MT. This belief 
is supported by the morphology: the needle-like shape closely resem- 
bles that known for a martensite colony 19 . Indeed, outside of the 
needle-like regions, all the blocks have the same SAED pattern and 
they apparently belong to one matrix grain and retain the original 
orientation. 

Another supporting evidence is the character of the boundary 
between the needle colonies and the matrix. These needles run in 
two different orientations in this original a-grain (Fig. 2e). For each 
needle, the boundaries are comprised of some straight lines and 
steps. The straight parts are edge-on {112} bcc twin boundaries, and 
on the needle (twin) side the {112} T plane (the green dotted line in 
Fig. 2a) is parallel to the boundary. But on the matrix side the {1 12} M 
plane (the black solid line in Fig. 2a) has a small tilt angle (cp) of —5° 
inclined with respect to the boundary. The {112} boundary itself is 
thus not a coherent twin boundary. These features are clearly seen in 
Fig. 3, from the high-resolution TEM image and the small angle 
between the {112} diffraction spots from the twin and the matrix 
in the corresponding SAED pattern in the inset. This character 
((p~5°) at the boundary is in fact expected from the compression- 
shuffle mechanism mediating the a^s transformation 9 under 
shock-loading (different from the hydrostatic pressure case 22 ), as 
schematically illustrated in Fig. 4a. A uniaxial loading compresses 
the atoms in the (110) bcc plane along the [001] bcc direction by 18.4% 
with respect to the uncompressed bcc state (top in Fig. 4a) to form a 
hexagon. Then every other (1 10)b cc plane shuffles by 3. 0 / (3^2) toward 
the [-110] direction (small arrows) to form an hep structure (bot- 
tom). The resulting cp between the resulting (l-100)h cp plane and the 
original (-1 12) bcc plane is 5.3°. We indeed observe cp = 5.2±0.7° stat- 
istically (for the statistical details, see Supplementary Table SI) 
between the (-112) bcc (as {112} M in the original a-phase in Fig. 2a) 
and the (l-100) hcp in the s-phase (as {112} T in the eventual a-phase 
in Fig. 2a). 

The mechanism of the s^a MT is illustrated in Fig. 4b. Shearing 
takes place along [11-20] hep (top row) or [-1-120] hep (bottom row) in 




Figure 3 | The special character of the needle-like colony's boundary 
revealed by HRTEM image. This particular example exhibits a —4.3° 
angle between the (-112) T and the (-112) M planes, which is also confirmed 
by the SAED pattern in the inset. 
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Figure 4 | Atomic mechanisms of the a-^>£-^>a phase transformation, (a) The a— »e phase transformation by compression-shuffle mechanism under 
shock-loading 9 , (b) The s^-a phase transformation by shear- shuffle mechanism. The (110) a and (000 l) e planes are in the plane of the paper. The red and 
yellow balls with bonds denote two atomic layers in or above the page, respectively. In (a), the bcc structure (top) is compressed along the [001]b cc 
direction by 18.4%, and every other ( 1 10) bcc plane shuffles by a 0 /(3j2) toward the [-1 10] bcc direction (small arrows) to form an hep structure (bottom). 
The [001] bcc compression without lattice expansion in the [-110] bcc direction results in a tilt angle of (p = 5.3°, between the resulting (l-100) bcp plane and 
the original (-112) bcc plane. In (b), shearing takes place either along [1 1-20] hcp (top row) or [-l-120] hcp (bottom row) in the (l-100) hcp plane, resulting in 
bcc lattices with a (-112) twinned relationship (compare iii and vi). Small shuffles (arrows) are required for various types of atoms (labeled 1-4 in the top 
middle panel). Executing the above shear-shuffle from the three crystallographically equivalent {l-100} hcp planes, three variants (V x through V 3 ) are 
possible for the eventual transformation product, resulting in the three sets of { 1 12} bcc twins that retain three-fold symmetry. 



the (0001) hcp plane. Meanwhile, the lattice parameter c 8 expands 
from 3.93 to 4.05 A while a s increases from 2.45 to 2.48 A, and the 
60° angles on the (0001) hcp plane change to the 70.56° on (112) bcc in 
the process. Near the shearing atomic row, atoms shuffle with various 
displacements: e.g., atom 1 through 4 in Fig. 4b(ii) by 0.832 A, 0.032 
A, 0.4258 A, and 0.454 A, respectively. Shearing in two opposite 
directions [Figs. 4b(i) vs. 4b(iv)] is equivalent geometrically but 
results in two mirror bcc lattices with regard to the (-112) bcc plane 
[compare Figs. 4b(iii) and 4b(vi)]. Executing the above shear-shuffle 
from the three crystallographically equivalent {l-100} hcp planes, 
three variants (V x through V 3 ) are possible for the eventual trans- 
formation product, resulting in the three sets of {112} bcc twins that 
retain three-fold symmetry. This explains how the s^a reverse 
transformation leads to the twinned structure inside the needle col- 
onies. What we actually observe under TEM thus matches exactly 
what one would expect to see in the microstructure from this s-^a 
MT mechanism. 

Discussion 

It should be noted that the shock pressure direction would be along a 
different orientation for each grain in our polycrystalline sample. In 
some original a-phase grains, the phase transition needles form pre- 
dominantly along only one {112} bcc plane (Supplementary Fig. S3), 
rather than along two crystallographically equivalent {1 12} bcc planes 
as in Fig. 2. Along the other {112} bcc planes, large Schmid factor 
activates plastic deformation (deformation twinning without the 
characters of the three-fold symmetry and the (p~5° at the twin 
boundary). However, all the characteristic features of the needle-like 
regions (transformed regions) in Supplementary Fig. S3 are the same 
with that in Fig. 2, so the same phase transition mechanism (Fig. 4) is 
activated irrespective of the varying grain orientations. 

Some researchers have predicted that a face-centered cubic (fee) y- 
Fe phase would occur under shock-loaded iron, based on molecular- 
dynamics simulation or numerical simulations 13 ' 21 . Here we have 
considered this possibility but excluded the involvement of inter- 
mediate y-Fe in the phase transition of shock-loaded iron. If the 
fee y-Fe was the intermediate phase instead of the hep s-Fe in the 



phase transition, the martensite a-Fe and the matrix y-Fe would have 
an orientation relationship of {lll} fcc //{011} bcc 23,24 . The {011} bcc 
planes of the martensite variants (lamellar twins) would present 
the angle of the {lll}f cc planes. The expected angle is 70.5°, which 
is not observed in our experiments. 

We now summarize the entire a-^E-^a process during shock 
treatment, using a schematic illustration in Fig. 5. In the original 
a-grain (green in Fig. 5a), the s-needles (blue) form along one or 
two {112} bcc planes during shock-loading (Fig. 5b), with a character- 
istic tilt angle at the colony boundaries, resulting from the compres- 
sion-shuffle mechanism known for the a^e MT. Afterwards, the 
needle-like s-phases transform back to the a-phases during unload- 
ing: three variants of the a-phase emerge inside the needles, inher- 
iting the threefold symmetry from the {l-100} hcp planes of the s 
structure. The shearing and shuffling mechanism involved in the 
e-^a transformation results in twinned structures, and the final 
configuration is therefore three {112}<111> twin sets (Fig. 5c) with 
an unusual 60° impinge angle. These morphological and orienta- 
tional features observed in the microstructure then serve as a tell-tale 
sign of the reversible MTs (a^e^a). 




Figure 5 | Schematic illustration of the a-^£^a phase transformations. 

(a) The original a-grain. (b) The a^>s phase transition under shock- 
loading, forming needle-like s regions, (c) The z-^ol reverse transition 
inside the s-needles upon unloading. Green and blue colors stand for the 
a-phase and the s-phase, respectively. Black lines are boundaries of the 
microstructural features. 



SCIENTIFIC REPORTS | 3 : 1086 | DOI: 1 0.1 038/srep01 086 



4 



Importantly, having identified clear microstructural indicators, 
the volume fraction transformed, F, can now be assessed not only 
quantitatively (the previously reported estimate under shock-loading 
had to fit an empirical function to the experimental Hugoniot 25 ), but 
also at different locations for a large sample after massive shock 
treatment. F can be estimated by directly measuring the area of 
needle-like regions on the EBSD micrographs acquired in scanning 
electron microscopy (SEM) at a number of locations along the radius 
in the whole upper surface of the shock-loaded hollow cylinder. As 
the distance inward from the loading surface increases, the shock- 
loaded pressures decreases, thus F decreases: for example, 54.6 ±3.2% 
(for the statistical details, see Supplementary Table S2) at 1.2 mm 
with -14.2 GPa (such as Fig. Id) and 27.1 ±2.8% (see Supple- 
mentary Table S3 for details) at 1.9 mm with —13.3 GPa (such as 
Fig. lc). Such a quantitative measure is very important for profiling 
and monitoring the shock response of a bulk material, and for val- 
idating models and predictions. 

In summary, in shock-compressed Fe we have detected and dec- 
iphered microstructural fingerprints that can be traced to ot-^s-^a 
transformation, using postmortem samples under TEM. While there 
are other in-situ measurements that can detect this phase transition 
under (local) shock-loading, our method provides more detailed 
information on the transformation mechanism. More importantly, 
the evidence we acquired is based on microstructures (characteristic 
morphologies). The signature features, i.e., the needle-type regions 
and their boundary characteristics as well as the angle between the 
twinned lamella inside, are uniquely revealing of the MTs involved. 
From them we can tell apart the transformed areas from those that 
have not experienced the MT. This has enabled a quantitative evalu- 
ation of the volume percentage of the transient s-phase formed at 
different locations and pressures under a complex shock-loading 
condition in a polycrystalline sample. In particular, our postmortem 
"microstructure" approach has advantages in removing many dif- 
ficulties and limitations associated with the use of relative intensities 
of in-situ X-ray diffraction peaks (e.g., signal to noise ratio, average 
over the probed thickness, limited sample shape/size/thickness and 
shock-loading conditions). The method of crystallographic analysis 
in the present work reproduced the formation process of the lamellar 
twins in the shocked iron. The same method may also be useful for 
inferring the origin of the Neumann bands which are mechanical 
twin lamellae formed under extraterrestrial collisions in bcc struc- 
tured kamacite phase in meteorites 26,27 . 

Methods 

Experimental setup. The experimental setup is schematically shown in 
Supplementary Fig. S4. Polycrystalline Fe with a purity of 99.90% was used to 
fabricate a 60 mm high hollow cylinder with a 42.5 mm outer radius (R) and a 
22.5 mm inner radius (r) (Supplementary Fig. S4a). The outside surface of the 
cylinder was covered with a 2.5 mm organic glass buffer layer followed by a 15 mm 
PETN explosive layer (as shown in Supplementary Fig. S4c). The entire explosive was 
detonated nearly at the same time by the detonating system illustrated in 
Supplementary Fig. S4c. The shock wave proceeds in the radial direction of the 
cylinder labeled by the purple arrows in Supplementary Fig. S4b. A VISAR equipment 
located in the center of the hollow cylinder (shown in Supplementary Fig. S4c) was 
used to record the inner free surface velocity (Supplementary Fig. S5). The structural 
evolution was investigated using the EBSD technique on a JEOL SEM 6700 with a field 
emission gun and TEM (JEM2010 and FEI Tecnai F30). 

Pressure calculations. The experimental setup cannot record the velocity for every 
point inside the sample. So the pressure history in the sample was calculated using a 
mixed phase model proposed by Hayes 28 and Andrews 29 based on a 2-D 
hydrodynamic code. Also, Boettger's model on metastable surfaces for forward and 
reverse phase transitions in pure iron and the parameters from his article 25 were used 
in the calculation. According to the calculation for every point along the radial 
direction in the sample, with increasing distance inward from the loading surface, the 
peak pressure decreases. At every point the explosively driven Shockwaves rise rapidly 
in pressure and then release at a slower timescale. The peak pressure of the location 
near the cylinder outside/shock loading surface is —16.4 GPa. To compare with the 
experimental observations, the calculated pressure histories of the points at 4.6, 1.9, 
and 1.2 mm inward from the loading surface are shown in Supplementary Fig. S6a. 
The peak pressures at the points 4.6, 1.9, and 1.2 mm are —12.7, —13.3 and 



~14.2 GPa, respectively. The inner free surface velocity curves from the calculation 
and experiment are shown in Supplementary Fig. S6b, which are in good agreement 
with each other. 
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